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Abstract: We investigate the mutiphoton process between different Bloch 
states in an amplitude modulated optical lattice. In the experiment, we 
perform the modulation with more than one frequency components, which 
includes a high degree of freedom and provides a flexible way to coherently 
control quantum states. Based on the study of single frequency modulation, 
we investigate the collaborative effect of different frequency components 
in two aspects. Through double frequency modulations, the spectrums 
of excitation rates for different lattice depths are measured. Moreover, 
interference between two separated excitation paths is shown, emphasizing 
the influence of modulation phases when two modulation frequencies 
are commensurate. Finally, we demonstrate the application of the double 
frequency modulation to design a large-momentum-transfer beam splitter. 
The beam splitter is easy in practice and would not introduce phase shift 
between two arms. 
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1. Introduction 

Ultracold atoms in periodically driven optical lattices, including shaken and amplitude modu¬ 
lated systems, have attracted much attention recent years, for they can bring out new interesting 
phenomena including realization of artificial gauge fields in different lattice geometries DHl] 
and the coherent control of atomic wavefunctions iHii. In a shaken one dimensional optical 
lattice, the ferromagnetic transition in trapped Bose gas has also been observed Q by coupling 
ground band s of the lattice system to the first excited band p. Specific to amplitude modu¬ 
lated lattice systems, although only bands with the same parity can be coupled, there are still 
researches on a wide range of problems, such as transfer of atoms from the ground band s to the 
second excited band d JS), detection of superfluid-Mott insulator transition © and study of the 
dynamical tunnelling of ultracold atoms with quantum chaos cniiii]. The technique can also 
be applied in areas including realization of a velocity filter IfT^ and detection of gravity 11311411 . 

Normally, studies of amplitude modulated lattices are based on single frequency modulations 
that only involve processes with one photon emission or absorption. During a polychromatic 
modulation, not only the amplitudes, but also the phases of different frequency components in 
the modulation can be controlled independently, providing a more flexible way to coherently 
manipulate quantum states. In this paper, we coherently transfer atoms from the ground state s 
to the high excited g band via a double frequency modulation. The peaks of transfer rate with 






Fig. 1. A sketch of lattice depth modulation in our system. The depth of lattice potential 
Vocos^{kix) is driven by a polychromatic modulation Vicos((n,f + <l>i). 


different lattice depths are measured experimentally, while influence of the modulation phase 
is demonstrated by performing an interference between two independent paths of excitations. 
These experiments completely investigated how frequencies and phases of different modulation 
frequency components would influence the excitation between Bloch states. Furthermore, we 
show an application of the double frequency modulation to build a large-momentum-transfer 
(LMT) beam splitter. Comparing with LMT beam splitters based on Bloch oscillation ISl fTSl or 
high order Bragg scattering ifT^ . our method is easy in practise and would not introduce phase 
shift between two separated atom clouds. 

The paper is organized as follows: In Sections |2]and[3the polychromatic modulation theory 
based on Floquet method and our experimental system are introduced respectively. In SectionH) 
we study the effects of the modulation phase for the single frequency modulation. In Section|5] 
collaborative effect of different frequency components is studied, demonstrate the effect of both 
modulation frequencies and phases. Section |6] presents a way to realize a LMT beam splitter. 
Discussion and conclusion are in Section [T] 

2. Theory for optical lattice with polychromatic modulation 

For an atom in an amplitude modulated lattice system along x axis, as schematically shown in 
Fig.m the time dependent Flamiltonian can be written as 

p2 

H{t) = ■^+Vocos^{kLx)+Y^ViCos{o)it + (l)i)cos^{kLx). (1) 

The first term in the right hand side is kinetic energy with M the atom’s mass and px its momen¬ 
tum along the x direction. The second term represents optical lattice without the modulation. 
Vo is the constant part of the lattice depth, and the wave vector is ki = 2njX with X the laser 
wavelength. The last term expresses the amplitude modulation with the modulation amplitude 
Vi, the frequency cu, and the phase (/),■ of each frequency components. 

Typically, time-periodic systems are described by Floquet’s theorem Oil . In our system, each 
of the modulation frequencies cu, gives a period 7] = InjcOi and the Hamiltonian has a period 
of H{t + T) — H{t), with T the lowest common multiple of 7;. By defining a time-evolution 
operator for one period as U[T], solutions to this problem must satisfy 

Ir,.a(f + F)) = f/(7’)|v/,,a(f)) (2) 

and the Floquet states \uq_a) are defined as with the quasi-momentum 

q, the band index a and quasi-energy «, which leads to 

{H[t) — idt)\uq^a) = Holuq^a) = £q,a\ucj,a) ■ 


( 3 ) 









For a special case of single frequency modulation, considering the time and coordinate peri¬ 
odicity of the Floquet states \uq,a), by a Fourier transformation we can study the problem with 
a set of basis |v/,„) = jn an extended Hilbert space ^ ® . Where M" is the 

Hilbert space and is the space of all functions with periodic T. By integration in one period, 
we can turn to a time-independent system. 


1 /■ 2 - 
1 Jo 


( 4 ) 


To get rid of the phase , we perform an unitary transformation [//„ 1^ = to the Hamil¬ 
tonian by A' = U^^AU. There are three kinds of terms in the matrix of A'. The diagonal 
terms / 2 M + mh(Oi are the energy of momentum states shifted by absorb¬ 
ing or emitting Floquet photons with energy Ticoi. The terms ^ = Vo/4 show stationary 

component of the lattice would coupling atoms with momentum difference 2 hki. In addition, 
mil “ terms show the time modulation of the lattice, which would coupling two 
momentum states separated by 2 hkL while absorbing or emitting one Floquet photon. 

When the modulation frequency is near-resonant with the energy difference between two 
specific bands and far detuned from the others, we can get an effective Hamiltonian by a rotating 
wave approximation. 


Hr 


Ea \ 

Vi.-'^>n*„p Ep-Tico J 


( 5 ) 


where Ea and Ep are the energy of Bloch states |a), |j 3 ) in the system without modulation. 
Coupling constant between two states is Q.ap = {a\cos^{k[jc)\P)/ 4 -. 

The extension to polychromatic driven is straightforward. For example, a double frequency 
modulation induced two-photon process between s-g band is described by an effective Hamil¬ 
tonian Hsg as; 


Hsg = 



( Es 


0 

0 

0 \ 




Ej —hcoi 

0 e‘1>^-V2£ldg 

e“^'Viadg 

0 




0 

Ed-hoh d'l’^Vikldg 

0 

e’^-V2£ldg 

.(6) 


0 


Eg-h{(0i+C02) 

0 

0 


0 


0 0 

Eg — 2hc0i 

0 



1 0 

0 

e-''^-V2El*g 0 

0 

Eg — 2h(l}2 y 


The effective Hamiltonian Hsg 

is constructed by means of nearly degenerate perturbation 


technique ifTSl . in which we include six nearly degenerate states considering four main pro¬ 
cesses in the excitation. The six states are \Es) the s band, \E^ — ticO\), \Eij — tioyi) the d band 
dressed by Floquet photon ( 0 i or 0)2 and \Eg — h{( 0 i + 0)2)), \Eg — 2 hcO\) and \Eg — 2/10)2) the 
g band dressed by two Floquet photons. Using this basis a general state (vi,v'2, V3, V4, V5, 
gives complex coefficient of the six dressed states. Population of Bloch states i is |vi p, while 
population on g band is |v4e4®i+®2)' _|_ yjgZirair _|_ given by coherent superposition of 

all g band states dressed with different Floquet photons. Solution of the model consists with the 
time dependent Schrodinger equation and the effective model provides us a better understand¬ 
ing of the mutiphoton process. However, in the calculation more states associated with higher 
order processes could be included to get a more accurate result, especially when the modulation 
amplitude is large. 
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Fig. 2. Left side is the calculated Floquet spectra of a single frequency driven system, 
with parameters Vq = 5Er, Vi = 0.5Er, TiC0\ = 5Er. In the figure the first seven bands are 
presented. The heavy lines depict states maximally overlapping with the i(blue), p (green) 
and d (red) Bloch bands respectively. Right side shows the details of two Floquet bands 
most overlapping with i and d bands. The two bands are separated by a band gap Ef at 
q = Q. 

3. Experimental system 

Our experiment begins with a quasi-pure condensate of typically 1.5 x 10^ Rb atoms in the 
jE = 2 ,ot/7 = 2) hyperfine ground state, produced in a combined potential of a single-beam 
optical dipole trap and a quadruople magnetic trap. The trapping frequencies are cp, =2n x 
28//z, (Oy = 2n X 6QHz, (O- = 2n x IQHz- The optical lattice is formed by a retro-reflected red 
detuned laser beam, with lattice constant a = Xj2 = 426nm focused to a waist of 110/im. 
Density of atoms in our system is less than 5 x lO'^cm^^, and the mean-field interaction can be 
omitted to capture the main physical mechanism in the excitation II19II20I . 

The lattice depth is calibrated by Kapitza-Dirac scattering and the modulation of lattice depth 
is controlled by an acousto-optic modulator(AOM). The modulation amplitudes Vi and phases 
0 , are generated from a signal generator and the intensity of lattice laser is monitored by a pho¬ 
todetector. Experimental results are absorption images taken after 2%ms time-of-flight (TOE). 
Occupation number at different momentum states \2lhkL) (I is integral momentum index) can 
be given from TOE images by «; = Ni/N, with Ni the atom number at momentum state \2lhki) 
and N the total atom number. The initial state for experiment is prepared non-adiabatically with 
numerically designed sequence of lattice pulses Il2ni22ll . Eor experimental convenience, the 
lattice pulses are carried out with the same depth as the constant part of the modulated lattice 
potential Vq. Typically, each of the pulses and the subsequent intervals are lasting for no more 
than 25/Xs, and the whole process can be finished within 60/is, thus the loading time is greatly 
reduced comparing with traditional adiabatic loading method. 

4. The initial phase effect in single frequency modulation 

Single frequency modulation can be seen as the basis of polychromatic driven. In this part, we 
present the preparation of a Eloquet state in the single frequency driven system, which is shown 
to be highly related to the modulation phase. 

Eollowing the discussion in Sec. |2] Eig. |2] depicts a typical quasi-energy spectrum of the 
single frequency driven system, which is obtained by direct diagonalization of A' at various 
quasi-momentum q. The same calculation also gives eigenvectors in the extended Hilbert space. 
The spectrum exhibits a complex structure as a result of the periodically repetition of high 
excited bands. 
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Fig. 3. Time evolution of n; measured from the experiments with initial modulation phase 
(al) ij) = —Till and (a2) ij) = ;r/2. Time averaged fraction (n;) are also shown for (hi) 
(j) = — ;r/2 and (b2) (j) = 7t/2 respectively, hq is shown with hlack dots comparing to the 
numerical simulation in solid lines. ni and are shown in average with red circles the 
corresponding numerical result is shown in dashed lines. Each point is averaged by three 
experiments and the error bars indicate the standard deviation. 


To connect our calculation with the experiments, we should project the Floquet states \uij^a) 
into momentum space. The occupation number at momentum states lllTiki) and its phase can 
be given by a summation of all the Fourier components with the same I as ci{t)\2lTikL) = 
where V/m = (v/m!«?,«) is coefficient of the Floquet state. T is related 
to the modulation phase and holding time f as T = j p is also useful to define the 

overlapping between the Floquet state and a Bloch state |«^) of the undriven potential 

as P = Jq I Property of the Floquet band is typically characterized by its most 

overlapping Bloch band ||23| . Without loss of generality, our experiments are restricted to quasi¬ 
momentum q = 0, and energy gap between Bloch bands a and j3 are written as ficOap ■ The 
technique can also be applied in systems with acceleration M, which brings out phenomena 
different from our study. 

When modulation phase 0i is given, a Floquet state can be projected into the momentum 
space, and such a state can be prepared by carrying out two lattice pulses with numerically 
designed pulse sequence 021II22I . In the fast loading process, lattice depth of the two pulses is 
kept constant. For a target state, the fidelity of a prepared state can be given numerically for 
different pulse durations and time intervals, and the optimized pulse sequence is obtained by 
finding the maximum loading fidelity. Throughout this method we can get a loading fidelity of 
more than 95% experimentally. 

With the initial modulation phase = —nj2, a pure Floquet state is loaded by lattice pulses 
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Fig. 4. (a) Two special cases in detecting the transfer population spectrum. In case 1, ab¬ 
sorption of photons with ®i (orange) or 0)2 (red) is resonant with d band. In case 2, two fre¬ 
quencies are equal, (b) For s-g coupling coi provides a two-photon process while Oh = 2a)i 
provides a one-photon process. Phases of two paths are controlled independently by mod¬ 
ulation phases of coi and o >2 ■ 


calculated with z = = —T/4. For the same state, we also present a modulation with 

= njl to show the influence of modulation phase. The population on momentum com¬ 
ponents \Qhki) and | -izlhki) for different initial phases = —7r/2 and = —njl are shown 
in Fig.[3al) and|3!a2) respectively. The parameters are Vq = S.OZsrjFi = 0.5£r and driven fre¬ 
quency is /toil = 5.QEr, red detuned from the band gap TicOsd = 5.2Er. Figure|3al) is a Floquet 
state, which shows a time period of T = 2n/cOi = 62.75/is. The time evolution in momen¬ 
tum space would be quite different when the modulation phase is changing by n, as shown in 
Fig.|3^a2). Besides the oscillation with a frequency near (Ui, we observe a slow growth of «i. 
Meanwhile, the amplitude of oscillations are also increasing following increase of the popula¬ 
tion on d band. 

The influence of modulation phase is shown more clearly when taking time average in one 
modulation period T. In Fig. [3bl) with = —7r/2 the Floquet state shows a constant pop¬ 
ulation in different momentum states by taking time average. While for = k/2. Fig. [3b2) 
shows a Rabi oscillation between s-d bands with Rabi frequency Ep/Ti, where Ep is the gap 
between two Floquet bands. The observation can also be explained by Eq. Q. The phase of 
coupling terms would change with 0i, thus the eigenstate is superposition of |s) and \d) with a 
relative phase determined by . When the modulation phase is changed, the loaded state is no 
longer an eignstate, and we can observe the Rabi oscillation. 

5. Excitation of ground state via a double frequency modulation 

The study of lattice modulation can be extended from single frequency to a more general form 
as described in Eq. (|6]l. Specific to two-photon excitation between s-g band, the time dependent 
lattice is described as 14(1) = Vq-\-V\ cos(a)if -F 0i) + V 2 Cos(a) 2 f + (j) 2 ), which includes seven 
parameters Vq; V'l , V 2 : , 0)2,0i, 02- During the experiments, we mainly focus on the influence 

of modulation frequencies and phases, while keeping other parameters constant. 

5.7. Spectrum of two-photon excitation 

The single frequency modulation shows that 0i is related to relative phase between s and d 
band components of the Eloquet state. However, with a nearly pure s band, the phase of single 
modulation can be neglected, only the phase difference between two modulations is important. 
Therefore, we chose (j)i = 7t while leaving 02 variable to control the relative phase between two 
modulations. 

Eor s-g band coupling through a two-photon absorption process, the sum of the two mod¬ 
ulation frequencies is chosen as Oi + 0)2 = (Osg- This process is well described by Eq. (|6l), in 













0.5 



20 25 30 35 40 



26 28 30 32 34 36 38 


aJi(kHz) 

Fig. 5. Spectrum for the population on ±ATikL states with increasing of modulation fre¬ 
quency (Bi. Population detected on ±4hki after a double frequency modulation for (a) 
Vb = 5£r(black) with Vi = lAEr, V 2 = 1.6Er t = 300^s, (b) Vq = 10£r(blue) with 
Pi = 2.SEr, V 2 = 2.2Er and t = 200ns, (c) Vo = \AEr{Kd) with Vi = V 2 = 2.5Er and 
t = ISO^i" are shown in rectangles with error bars. Solid lines are corresponding numerical 
simulation. 


which we have considered different pumping paths, as schematically shown in Fig.|4ja). There 
are two cases which would benefit the excitation process. 

Case 1: Resonant two-photon process. When 0}\ — CO^d or (02 = COsd^ atoms are transferred 
from |s) to |g) with the assistance of d band as an intermediate band. 

Case 2: Equal frequency two-photon process. When ( 0 i=O >2 = (Osgl^., two modulations with 
the same frequency can be added together, and the coupling strength of the process is doubled. 

With resonance condition Oi = (Osd in case 1, when two modulation amplitudes are cho¬ 
sen as Vi = V2^dg the transfer rate would show a maximum resonant peak. And we keep 
this modulation amplitudes while sweeping (Oi . Modulation phase is chosen from numerical 
simulation to get a maximum transfer. 

In the experiments, we sweep the frequency coi for different lattice depths Vq = 5Er, IQEr and 
lAEr- Population on momentum states ±Ahki measured from the experiment are shown with 
rectangles in Fig.|5] comparing with the theoretical calculation shown in solid curves. Within 
the lattice depth we considered, g band is greatly concentrated on | ± ATiki) momentum states, 
thus n ±2 can reflect transfer rate to g band. 

Figure 13 a) shows the case of Vq = SF, - For the lattice depth we have Ei^d/Eidg = 1.11, 
correspondingly the modulation amplitudes are chosen as Vi = lAEr, V 2 = 1.6Er, and t = 
300/is. The holding time t may be chosen shorter if the maximum of numerical simulation is 
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Fig. 6. The excited population on g band shows the interference between two paths, (a) 
Population transferred to n ±2 is shown in black dots with error bars. The dashed line shows 
theoretical simulation for comparison. (bl)-(b4) Vi for different phases. 


reached at an earlier time. In the figure, two peaks appear at Oi = (Osd ft )2 = (Osd which 
follows case 1 we have discussed. And there the central peak at frequency (£)i= (O2 following 
case 2 is much lower than two peaks for case 1. 

Figure|5lb) shows the spectrum with Vq = lOEr- With the increasing of Vq, the energy differ¬ 
ence (Osd is getting closer to cOdg, and three peaks are overlapping. Comparing with Vq = 5Er 
the central peak is much higher, because the process of case 2 is also near resonance with d 
band. 

For Vo = l4Er only one peak would be measured in the spectrum as shown in Fig. He), 
which means easel and case2 are fulfilled simultaneously. Under this condition, the coupling 
between s and g band is also greatly enhanced. Modulation amplitudes are Vi =V2 = 2.5Er, for 
the two modulation frequencies are the same at easel and can’t be distinguished. 

The experimentally detected peaks are governed by two cases, which are within the descrip¬ 
tion of Eq. (|6]l. Thus in the numerical simulation we can neglect higher order processes of 
emission and absorption of Floquet photons. The discrepancy between experimental result and 
the theoretical simulation is probably due to the influence of interaction and initial momentum 
distribution of the condensate. These effects would destroy coherency during the modulation, 
and the measured population of excited state would be lower than the maximum value in the 
theoretical simulation. 

5.2. The role of modulation phases 

When <Oi and (O2 are incommensurate, the influence of relative phase is not prominent because 
there is only one path for the pumping process and no states could interfere with each other. 
However, the effect of modulation phases would be more pronounced when two modulation 
frequencies are commensurate. 

As shown in Fig. HJb), when 2(Oi = © 2 , relative phase of modulations can be shown by 
performing a one-photon process simultaneously with a two-photon process. We choose the 
frequency coi = ©vg/2 at central peak of Vq = IQEr lattice and the second modulation is per¬ 
formed with ©2 = (Osg. Modulation amplitudes are Vi = V 2 = 2.5Er and t = SOOps. Similar 
to the process described by Eq. ©, in this problem we consider the interference between two 
states \Eg — 2h(o{) and \Eg — hcof). During the two-photon process with 0i = n, phase of state 
\Eg — 2h(Oi) remains zero, while 02 determine phase of state \Eg — hcof) as 02 — 7t. Relative 
phase of two modulations is defined as (02 — ;r) — ^(0i — ;r) = 02 — K, and the interference 
can be changed from constructive to destructive with different 02 . 
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Fig. 7. A LMT beam splitter with a separation of Hhki. (a) TOF image of the LMT beam 
splitter, (b) Experimentally measured population of atom on momentum states | ± 6hki} 
are shown in black dots with error bars. 

Figure |6ja) depicts the population transferred to g band with the phase of one-photon pro¬ 
cess changing by Itt. Figure|6lb) shows how the depth of lattice is varying with time for four 
different modulation phases. With 02 = TT, the population transferred to |g) through one-photon 
process and through two-photon process are in phase and the transfer rate is enhanced. When 
02 is increasing, relative phase of two processes are deviating from zero, and the population on 
|g) would decrease. In the case of 02 = 0 the two process are out of phase, and only few atoms 
can be transferred to g band. Further increasing of 02 would increase the measured population, 
and when 02 is changed by 27r the population on g band reaches the maximum value again. 

The single atom Hamiltonian in Eq. ([TJ can well explain the excitation in our system. How¬ 
ever, for the SQQjis modulation, decoherence from mean-field interaction and the momentum 
distribution becomes significant, thus it is necessary to carry out a simulation based on the 
time-dependent Gross-Pitaevskii equation. 



(7) 


where g is the parameter of interaction. In the simulation, we consider both the distribution of 
initial momentum and the mean-field interaction which would deduce the measured excitation 
rate. Dashed line in Figlbja) gives the result of simulation which fits well with the experiment. 
Similarly, the relative phase can also be changed by 0i of the two-photon process, which gives 
a period of k. 

The collaborative effect of different frequency components is studied in two aspects, demon¬ 
strate the effect of both modulation frequencies and phases of different frequency components. 
The usefulness of this study is not limited to double frequency modulation, within these two 
aspects, collaborate effect of more frequency components can also be well understood. 

6. Application in LMT beam splitter 

The double frequency modulation can also be applied to realize a LMT beam splitter by pump¬ 
ing atoms to higher momentum states, which is useful in experiments of atomic interferometry 
ll24l . According to the resonant condition we have discussed, a preferred choice is Vq = 14-Er 
with the modulation frequency COi = COsd = COdg- The other modulation frequency is chosen as 
0)2 = (Ogi (with i the 6'^ excited Bloch band) to get a distribution at ±6hki. 

We begin with a condensate, the lattice is suddenly turn on, modulated with Vi =¥2 = 2.5Er 
and the preferable phases are found numerically. A TOF image of experimental result is shown 
in Fig.|2ta). Figure|7jb) shows the population of atoms on momentum states | ± b/iki) where we 




have subtracted the thermal gas. Near 80% of the atoms are coherently transferred into ±6tikL 
momentum states within 160/is, and the maximum lattice depth needed is below Vo +Vi + V 2 = 
24£’r- Comparing with a LMT beam splitter based on high order Bragg scattering ifThl . our 
method needs a much lower lattice depth and is easy in practice. Furthermore, the process is 
symmetric for both sides, and would not introduce phase shift between two separated atom 
clouds. 

A momentum splitting of 12hki is not the limit of the amplitude modulation. More frequency 
components or another subsequent modulation can be introduced to reach a larger momentum 
splitting. For example, after the double frequency modulation, by preforming another single 
frequency modulation resonance with the energy difference between | ± 6hkL) and | ± ^Tiki) 
momentum states, the atoms can be transferred to | ± Shki) coherently. 

7. Discussion and conclusion 

Following Floquet-Bloch theory we have presented an experimental preparation of a Floquet 
state and study its property. The non-adiabatic loading method would take much less time and 
greatly reduce the heating problem. The loading method could also be extended to systems with 
a shaken lattice Q or a combined modulation |l 6 ]| 23 , which provides a several millisecond 
longer lifetime for condensate in the experimental study on areas including the detection of 
Floquet topological states 017ll251l26ll . 

In conclusion, based on the study of single frequency modulation, we investigate double 
frequency modulation in detail. Experimental observations show that different modulation fre¬ 
quency components would influence each other in two ways. When two frequencies are reso¬ 
nant with subsequent excitation processes, the modulation can induce a resonant two photon 
process which can effectively transfer atoms to higher excited bands. With specific modulation 
frequencies, interference between a one-photon path and a two-photon path is observed, re¬ 
vealing influence of modulation phases when two frequencies are commensurate. The quantum 
interference can be used to enhance the excitation rate or destruct unwanted excitations. Using 
the technique of double frequency modulation, we also demonstrate an efficient way to realize a 
LMT beam splitter with low lattice depth. Our study provides a more flexible way to coherently 
control quantum states through an optical lattice. 
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